A pQCD-based description of heavy and light flavor jet quenching 
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We present a successful description of the medium modification of light and heavy flavor jets 
within a perturbative QCD (pQCD) based approach. Only the couplings involving hard partons are 
assumed to be weak. The effect of the medium on a hard parton, per unit time, is encoded in terms 
of three non-perturbative, related transport coefficients which describe the transverse momentum 
squared gained, the elastic energy loss and diffusion in elastic energy transfer. Scaling the transport 
coefficients with the temperature of the medium, we achieve a good description of the centrality 
dependence of the suppression and the azimuthal anisotropy of leading hadrons. Imposing additional 
constraints based on leading order (LO) Hard Thermal Loop (HTL) effective theory leads to a 
worsening of the fit, implying the necessity of computing transport coefficients beyond LO-HTL. 
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Experimental results from the Rclativistic Heavy Ion 
Collider (RHIC) have established that a new kind of 
hot and dense partonic matter has been created in cen- 
tral Au+Au collisions [l|. High transverse momentum 
(pr) partons created in the initial hard collisions were 
predicted to provide a reliable probe of this highly ex- 
cited matter [2j. These partons lose energy in the dense 
medium leading to a depleted yield of high pt hadrons 
compared to that in binary scaled p+p collisions [3] . Due 
to the large energy scale involved, these hard partons 
were expected to couple weakly with the medium allow- 
ing for the use of perturbative QCD (pQCD) to describe 
their propagation through the dense matter. Experimen- 
tally, the basic picture of parton energy loss has been 
confirmed by the observation of significant suppression of 
the high px yield for both light 0] and heavy flavors Q ■ 
There now exist sophisticated (and successful) calcula- 
tions of light flavor suppression [H| . Not only have these 
accounted for the centrality dependence and azimuthal 
anisotropy of the single inclusive suppression but also for 
photon and hadron triggered correlations on both the 
near [|| and away side of the trigger hadron Q. 

The success of pQCD based calculations of heavy flavor 
modification, however, has been less than satisfactory. In 
a prior effort, Armesto et al. @ were able to describe both 
light and heavy flavor suppression in central collisions 
(and the azimuthal anisotropy) including only radiative 
energy loss. However, they required a time averaged jet 
transport parameter q = d(Ap±) 2 /dt ~ 14 GeV 2 /fm for 
a gluon jet (p± is the momentum trasverse to the jet axis 
and t is the time spent in the medium), yielding a p± 
comparable to the energy of the parent jet and at least 
a factor of five larger than estimates from LO-HTL us- 
ing a medium with an identical temperature profile. In a 
recent analysis by Wicks et al. 0, the authors incorpo- 
rated both radiative and elastic energy loss. The clastic 
energy loss coefficient e — dE/dt was calculated in LO- 
HTL but the radiative energy loss was calculated in a 



different medium of static scattering centers, estimating 
the density of scatterers and the transverse momentum 
per scattering with the help of LO-HTL calculations con- 
strained by entropy considerations. These authors were 
unable to quantitatively describe the suppression of non- 
photonic electrons in central collisions at RHIC. This has 
led the authors of Ref. [l(| to question the ability of any 
pQCD based approach to consistently describe high pt 
suppression, leading to the speculation that hard partons 
with energies in the tens of GeV may be strongly inter- 
acting with the produced dense medium [ll[. This has 
cast a pall of doubt on the entire program of pQCD based 
jet tomography of dense strongly interacting media. 

In this Letter, we provide a satisfactory explanation 
of both light and heavy flavor suppression (including its 
centrality dependence and azimuthal anisotropy) within 
a single pQCD based approach. Taken in combination 
with the results of Refs. 043; this effort closes the gap in 
the suite of hard probe observables at RHIC that can be 
successfully described by pQCD. We base our approach 
on the factorization paradigm inherent to the higher twist 
scheme, where a hard jet is weakly coupled with the gluon 
distribution of the medium, which itself may be strongly 
or weakly coupled. The effect of the medium is encoded 
in terms of three non-perturbative transport coefficients: 
q, e, and the diffusion in elastic energy transfer per unit 
time e 2 = d(AE) 2 /dt. 

The primary difference with Ref. 0] is that we do not 
insist on estimating the transport coefficients in LO-HTL 
but instead treat them as parameters of the model. As- 
suming near on-shell propagation for the hard partons, 
we relate e and 82 to the loss in longitudinal momentum 
per unit time (dp z /dt) and diffusion per unit time in lon- 
gitudinal momentum transfer (d(Ap z ) 2 /dt). Assuming a 
medium close to local thermal equilibrium and the ap- 
plicability of the fluctuation-dissipation theorem [l2| , we 
set d(Ap ± ) 2 /dt ~ 2d(A Pz ) 2 /dt ~ (4T/\v\)dp z /dt, with v 
the velocity of the jet parton. The only input parame- 
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tcr is chosen as q (assuming the same for both light and 
heavy quarks) which is assumed to scale with temper- 
ature (q = CT 3 ). This form is similar to the results 
obtained for a hard jet weakly coupled to a strongly 
coupled medium [l3| and also suggested by higher or- 
der calculations of q in HTL effective theory 14| . With 
such setting, we achieve a satisfactory description of the 
medium modification of both light and heavy flavor jets 
and their elliptic flows. For completeness, we also per- 
form a fit with the form of the coefficients taken from 
LO-HTL theory (with the in-medium a s as the fit pa- 
rameter) . While the transport coefficients obtained from 
the fit turn out to be comparable, a worsening of the fit is 
obtained, in agreement with Ref. [{J. Thus, while the ap- 
plicability of LO-HTL effective theory to the calculation 
of medium properties may be questioned, the assumption 
that a hard jet is weakly coupled with the medium, al- 
lowing for a pQCD based description of light and heavy 
flavor suppression, is further verified and established. 

High momentum partons produced in hard colli- 
sions, tend to be rather virtual. In vacuum, vir- 
tually (Q 2 ) is lost by subsequent emissions. The 
probability for these emissions can be calculated us- 
ing pQCD as long as the virtuality at a given emis- 
sion is large enough. Below a certain Q 2 , one will 
need to use an experimentally fitted fragmentation func- 
tion (FF) to calculate the inclusive distribution of sin- 
gle hadrons. The change of this FF due to prior higher 
virtuality (transverse momentum) emissions can be cal- 
culated using the Dokshitzcr-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution equations [l5[. These require 
the measured FF at some lower scale /i 2 as input and 
compute the modification to this single hadron distribu- 
tion due to multiple emissions from /i 2 up to Q 2 . 

Hard virtual partons traversing a dense medium will 
scatter off the constituents of the medium in addition 
to emitting lower virtuality partons; this changes the 
momentum distribution of the propagating partons be- 
tween emissions. A parton with light-cone momentum q~ 
(q~=q — q z , for a parton traveling in the — z direction) 
much larger than its virtuality q~ 3> Q, can be effec- 
tively described with a length {Lr ) dependent three di- 
mensional distribution <f>{q~ , q±, L~). The change in this 
distribution due to multiple scattering per unit len gth , 
up to second order in gradients, can be expressed as [161 ] 
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where, q\ c , e\ c and e2i c are non-perturbative transport co- 
efficients which encode properties of the medium. These 
can be expressed in terms of two gluon operators, e.g., 



q lc =[A7r 2 a s C R /(N 2 - l)]jdy-(F+^y-)F+(0)), (2) 
e 21c =[4ir 2 a s C R /(N 2 - l)]fdy-(F+-(y-)F + -(0)). 



These operator products can be factorized from the hard 
process and computed in any given model of the medium. 
The angled brackets () indicate an expectation of the 
operator product in an arbitrary ensemble. We do not 
present the operator expression for the elastic loss e\ c as 
we will always relate it to the expression for fluctuations 
in the elastic loss, e2i c , via the fluctuation-dissipation the- 
orem. These light-cone transport coefficients may be re- 
lated to the Cartesian coefficients as 

q = (1 + \v\)q\ c ,e = \v\e\ c , e 2 = u 2 /(l + \v\)e 2 i c ■ (3) 

The multiple scattering and emissions from a hard vir- 
tual parton and its effect on the final single hadron distri- 
bution can be calculated as long as the virtualities, both 
on entry and exit from the medium, are large enough for 
pQCD to be applicable. For the case of only transverse 
scattering, when Q 2 3> fc 2 ^, the soft scale of the medium, 
one recovers a DGLAP like picture with ordered emis- 
sions given by a splitting function which is modified by 
multiple transverse scattering in the medium [13]. The 
medium modified fragmentation function (MMFF), for 
a hadron with momentum fraction z of the jet, now de- 
pends on the location of entry and exit (Q) from the 
medium; its change with virtuality is given as, 
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In the equation above, P%^j is the in-medium splitting 
function (IMSF) at location ( [13] , 
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where T£ is the formation time of the emitted gluon with 
forward momentum fraction y and transverse momen- 
tum l± = Q: Tf = 2q~y(l — y)/l\- The full evolution 
equations for MMFFs will also include a pure vacuum 
contribution, which is implicitly included. By taking 
mass effects into account, the above results have also 
been extended to heavy quarks [3] , which suffer less ra- 
diative energy loss due to the dead cone effect. Here 
r f = 2g-y(l - y)/[l\ + (1 — y) 2 M 2 }, with M Q the heavy 
quark mass. The IMSF also contains a multiplicative 
factor Q 8 /[Q 2 + (l-y) 2 A/ 2 ] 4 . 

In this calculation, the masses of the heavy quarks 
and heavy mesons are taken as: M c = 1.3 GeV, M& = 
4.2 GeV, M D = 1.9 GeV and M B = 5.3 GeV. Our in- 
put to the evolution equations is picked using restric- 
tions based on formation length: given a mean length 
traversed of {L~), partons with formation lengths (r = 
q~ /Q 2 ) much larger than this will not be modified by the 
medium, as a result, we take as input the vacuum frag- 
mentation function at /iq = q~ / L~ . In cases where /Kq 
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falls below 1 GcV 2 , /ig is set equal to 1 GeV 2 . These frag- 
mentation functions arc then evolved up to the hard scale 
using a sum of the medium modified evolution equations 
[Eq. and the standard vacuum evolution equations. 
The modification due to elastic energy loss is incorpo- 
rated by shifting the fraction z, 

D'(z)= { dAzP(Az)D(z/(l-Az))/(l-Az), (6) 



where P(Az) is a Gaussian distribution with a mean and 
variance determined by e and e-2 respectively. 

Final hadron spectra are obtained as the convolution, 



d(Th = ^2 fa/A ® fb/B ® d(T ab ^ jd <g) D 



(7) 
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where f( a ,b)/(A,B) is the nuclear parton distribution func- 
tion (PDF), da is the partonic cross section, and -D/i/j 
is the FF after radiative and collisional energy losses. 
For both light and heavy sectors, PDFs are taken from 
CTEQ5 1 1 Oil (shadowing corrections are obtained from 
EKS98 |20|). The partonic cross section da is evalu- 
ated at leading order, with a if -factor accounting for 
next-to- leading order effects (1.7 and 2.4 for light and 
heavy sectors, repectively) . The renormalization and fac- 
torization scales are set to be the transverse energy of 
hadrons = (p£ + Ml) 1 ' 2 . The vacuum FFs for light 
quarks and gluons arc taken from KKP |2l| and c— >D, 
b^B FFs are obtained from PYTHIA6.4j22j]. We in- 
troduce multiplicative semi-leptonic decay functions of 
heavy mesons {fe/n) to obtain the non-photonic electron 
spectrum. They are obtained by fitting to D decay from 
BABAR Q and B decay from CLEO 0. The nor- 
malizations are determined from their branc hing ratios: 
BR(D-^e) w 10% and BR(B— >e) = 10.36% jHf. With 
these, we obtain a reasonable description of the baseline 
experimental measurements of 7r° spectra pij and non- 
photonic electron production [26[ in p+p collisions. 

Given the dependence of q, e and £2 on temperature 
T, a midrapidity thermal space-time profile T(tj_,t) is 
required. In this application, the initial transverse spatial 
profile for the entropy density s ~ T 3 is chosen to be 
proportional to the participant density in the collision of 
two nuclei with a Woods-Saxon nuclear density function 
(parameters are taken from Ref. [13] )• The medium is 
assumed to thermalize at proper time tq = 0.6 fm/c, 
with a temperature To = 400 MeV at the hottest point in 
central collisions. The temperature diminishes with time 
as r -1 / 3 due to one dimensional Bjorken expansion. The 
spatial distribution of the initial hard jets is determined 
by the binary nucleon collision density. 

We compare with the measured nuclear modification 
factor (Raa) as a function of pr, defined as the ratio of 
the yield in A + A to the binary scaled yield in p + p, 
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1 dN AA /d 2 p T dy 
JVcoii dNPP/d 2 PT dy ' 
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FIG. 1: (Color online) Raa (solid lines) for tt° in central (a) 
and mid-central (b) and non-photonic electrons in central (c) 
Au+Au collisions at RHIC, with q oc T 3 (and corresponding e 
and £2) . Dashed and dot-dashed lines represent the radiative 
and collisional contributions, respectively. 
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FIG. 2: (Color online) The predicted minimum bias V2 for 
7t° and heavy flavor electrons in Au+Au collisions at RHIC, 
using the same parameters as determined in Fig. [T] 



where A^ co ii is the average number of binary collisions be- 
tween the nucleons from two nuclei in the range of impact 
parameters chosen. With a single qo w 1.3GeV 2 /fm (and 
associated values of e, 62) for a quark jet at T = 400 MeV, 
a satisfactory description of Raa for 7r° in central and 
semi-central collisions and the supression of heavy-flavor 
electrons in central collisions is achieved [see Fig. [T] . It is 
clear that the inclusion of elastic energy loss [28( is nec- 
essary to explain the suppression of both light and heavy 
flavor suppressions. Using these values we can predict 
the impact parameter integrated azimuthal anisotropy 
[minimum bias V2 , (Fig. 12}]- A non-Glauber-bascd initial 
profile will affect this final predicted V2, as will a more 
realistic calculation in a three dimensional hydrodynamic 
simulation. 

We close with a discussion of a weakly coupled medium 
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HTL to compute jet transport coefficients. 

We thank Ulrich Heinz and Berndt Midler for helpful 
discussions. This work was supported in part by U.S. 
D. O. E. under grant Nos. DE-FG02-01ER41190 and 
DE-FG02-05ER41367. 
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FIG. 3: (Color online) Same as Fig. [TJ but with q, e and e.2 
calculated in LO-HTL. 



where the transport coefficients are estimated using LO- 
HTL theory. This gives, for a light quark or gluon, 



<Zhtl = C R a s m 2 D T In [AET/m 2 D ] , 



(9) 



where 



= Ana s (l + N f /6)T 2 is the Debye mass 
(Nf = 3). LO-HTL expressions for the heavy quark e 
are obtained from Ref. [29j . The in-medium coupling 
a s is adjusted to describe the central collision data for 
light flavors; this yields qo ~ 1.3 GeV 2 /fm for a quark 
jet with E = 20 GeV. Using these we predict the results 
for non-central collisions and heavy quark energy loss. 
In agreement with Ref. Q, we find (see Fig. [3]) a con- 
siderable worsening of the comparison with the heavy 
flavor data (x 2 /d.o.f. = 110/20 ps 5.5, compared with 
X 2 /d.o.f. = 24/20 ps 1.2 in Fig. Q]). 

In conclusion, we have carried out a global fit of high 
Pt single inclusive observables for both light and heavy 
flavors at RHIC, where only the coupling of the jet with 
the medium is assumed to be weak. The properties of the 
medium are encoded by non-perturbative transport coef- 
ficients q, e and e.2, which are related by the fluctuation 
dissipation theorem for an isotropic thermal medium. 
With only one input parameter (chosen as q) scaling with 
T 3 , we achieve a good agreement with the experimental 
measurements for both the suppression and the elliptic 
flow. This supports the assertion that perturbativc QCD 
and weak coupling approaches can be applied to the de- 
scription of jet modification in a dense medium, even if 
the medium may itself not be weakly coupled. Insisting 
on estimating all the transport coefficients solely within 
LO-HTL theory leads to a failure to explain the heavy 
flavor suppression, implying the need to go beyond LO- 
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